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Abstract: In order to cope with drought and water shortages, the working people in the arid areas of 
Northwest China have developed a drought-resistant planting method, namely, gravel-sand mulching, 
after long-term agricultural practices. To understand the effects of gravel-sand mulching on soil water 
evaporation, we selected Baifeng peach (Amygdalus persica L.) orchards in Northwest China as the 
experimental field in 2021. Based on continuously collected soil water stable isotopes data, we evaluated 
the soil evaporation loss rate in a gravel-sand mulching environment using the line-conditioned excess 
(lc-excess) coupled Rayleigh fractionation model and Craig-Gordon model. The results show that the 
average soil water content in the plots with gravel-sand mulching is 1.86% higher than that without 
gravel-sand mulching. The monthly variation of the soil water content is smaller in the plots with 
gravel-sand mulching than that without gravel-sand mulching. Moreover, the average lc-excess value in the 
plots without gravel-sand mulching is smaller. In addition, the soil evaporation loss rate in the plots with 
gravel-sand mulching is lower than that in the plots without gravel-sand mulching. The lc-excess value was 
negative for both the plots with and without gravel-sand mulching, and it has good correlation with 
relative humidity, average temperature, input water content, and soil water content. The effect of 


gravel-sand mulching on soil evaporation is most prominent in August. Compared with the evaporation 
data of similar environments in the literature, the lc-excess coupled Rayleigh fractionation model is better. 
Stable isotopes evidence shows that gravel-sand mulching can effectively reduce soil water evaporation, 
which provides a theoretical basis for agricultural water management and optimization of water-saving 
methods in arid areas. 
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1 Introduction 


China's total freshwater resources account for 6% of the world's total, and its per capita annual 
water resources account for only 28% of the world's total. The amount of water resources per 
hectare of arable land in China is 2.1x10* m°, accounting for 50% of the world average (Kang, 
2019; Yuan et al., 2019). The shortage of water resources has become a severe challenge to 
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China's agricultural development. In Northwest China, gravel-sand mulching emerged as a 
characteristic agricultural technology in agricultural activities. It is a drought-resistant planting 
method developed by people working in the arid environments after long-term farming practices 
adapted to drought and water shortage. Gravel-sand mulching refers to covering the soil surface 
with gravel or sand, thereby inhibiting soil water evaporation (Zhao et al., 2017) and increasing 
soil water storage (Poesen et al., 1997). Evaporation is the main way of soil water consumption 
and the main output of the regional water cycle (Che et al., 2020). It is also an important factor 
causing soil water loss and drought. Therefore, reasonable reduction of soil water evaporation can 
effectively increase the total crop water supply and comprehensively reduce the impact of water 
storage on crop growth. 

The traditional evaporation assessment mostly adopts micro lysimeter, soil heat pulse, chamber, 
micro-Bowen ratio energy balance analysis, and eddy covariance (Kool et al., 2014). However, 
these methods are time-consuming and expensive when obtaining the long-term evaporation data. 
Although modeling is useful for estimating the soil evaporation loss rate over a relatively long 
period, few models can estimate evaporation under the canopy; even if they can, they tend to 
overestimate the amount of soil evaporation loss rate, which needs to be verified by long-term 
field monitoring (Xiang et al., 2021). 

The stable isotopes ratio of water (i.e., 7H/'H and !*O/'°O) is a good tracer for studying water 
movement and flux (Kool et al., 2014; Sprenger et al., 2017). Evaporation of the soil surface is 
the main cause of soil water stable isotopes enrichment, and the degree of evaporation depends on 
the time of surface exposure as well as other soil and climatic conditions (Hasselquist et al., 2018; 
Ducklert et al., 2019). To date, a lot of study has been done to estimate the soil evaporation loss 
rate under experimental and field conditions through the isotopic composition of water in shallow 
soils (Barnes and Allison, 1984; Hsieh et al., 1998; Sprenger et al., 2017; Al-Ogaili et al., 2020; 
Mahindawansha et al., 2020). As a supplement to the traditional method, the main advantage of 
the method based on stable isotope is that the estimates correspond to a relatively long period 
(such as days and months) before soil sampling, and it does not require continuous field 
observations (Barnes and Allison, 1984). In the past, the classical Craig-Gordon model was 
usually used to evaluate the evaporation loss of soil or lake by the stable isotope method (Craig 
and Gordon, 1965; Gibson et al., 2016; Sprenger et al., 2017; Che et al., 2020; Mahindawansha et 
al., 2020). Utilizing the measured values of 5!8O and 67H in water, this model can be used to 
calculate the soil evaporation loss rate. However, the soil evaporation loss rate calculated based 
on ô!ŝO and 87H is usually different in practical applications (Sprenger et al., 2017; Hu et al., 
2018; Mahindawansha et al., 2020). To make the calculation results more accurate, we used the 
dual-isotope systems method (line-conditioned excess (Ic-excess) coupled Rayleigh fractionation 
model) and the traditional single-isotope systems method (Craig-Gordon model) to quantify the 
soil evaporation loss rate in this study. 

Baifeng peach (Amygdalus persica L.) is widely planted in Gansu Province of China, and the 
economic benefits are considerable. In the relatively harsh natural environment, the healthy and 
sustainable development of economic forests is particularly important. Therefore, the evaluation 
of soil evaporation loss rate in the gravel-sand mulching environment of Baifeng peach orchards 
can provide data support for agricultural production activities, which is of great significance for 
the efficient use of limited water resources to maintain the sustainable and healthy development 
of fruit trees and support water management in orchards. In order to determine the impact of 
gravel-sand mulching on soil evaporation, we monitored the changes of soil water content and 
soil water stable isotope values in the plots with and without gravel-sand mulching during the 
peach tree growing season (from April to October) in 2021. The objectives of this study were to: 
(1) estimate the effect of gravel-sand mulching on the soil water content; (2) assess the soil 
evaporation loss rate in the gravel-sand mulching environment; and (3) evaluate the reliability of 
the soil evaporation loss rate calculations. 
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2 Materials and methods 


2.1 Study area 


The experimental field is located in the western Loess Plateau of China (36°08'N, 103°40’E) and 
belongs to the arid and semi-arid climate zone (Dang, 2004). The average annual precipitation is 
312.9 mm, the annual average temperature is 7.4°C, the potential evaporation is 1676 mm, the 
sunshine duration is 2447 h, and the frost-free period is 180 d (Chen et al., 2013). The natural soil 
is mainly dark gray calcareous soil and typical gray calcareous soil, which is barren and poor in 
soil and water conservation. The arable soil is mainly anthropogenic, which is the result of 
long-term cultivation, irrigation, and fertilization. The experimental field covers an area of 666.7 
m2, with irrigation canals in the north and east, other farmlands in the south, roads in the west, 
and there are two ridges in the field (Fig. 1). 
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Fig. 1 Location of the experimental field in the Loess Plateau (a), overview of the experimental field (b), and 
photographs of the plots with gravel-sand mulching (c) and without gravel-sand mulching (d). DEM, digital 
elevation model; A, C, and E are the plots with gravel-sand mulching; B, D, and F are the plots without 
gravel-sand mulching. 


2.2 Sampling and measurement of soil 


We selected six well-growing peach trees (Fig. 1b) in the experimental field, distributed on both 
sides of the two ridges; all the trees are Baifeng peach. In Figure 1b, trees marked with A, C, and 
E are covered with 5.0-cm thick sand and gravel after soil sampling in May 2021, and trees 
marked with B, D, and F are the plots without gravel-sand mulching. The radius of coverage area 
is 0.8 m, and the particle size of sand and gravel is less than 5.0 mm. The coverage measures used 
in this study are based on the experiences of local farmers. Soil samples were collected at a 
distance of 0.8 m from the peach tree. After manually excavating the soil profile at 0.5 m, a soil 
drill with a depth of 1.0 m was used for deep sampling. The sampling interval is 0.1 m, and the 
total sampling depth is 1.5 m. We collected 4 parallel samples from each layer of soil, and a total 
of 2688 samples were collected. A part of the soil samples was placed into aluminum boxes to 
determine the soil water content. The other part was sealed in a glass bottle and then frozen and 
stored until experimental stable isotope analysis was carried out. 

All the collected samples were analyzed in the stable isotopes laboratory of the College of 
Geography and Environmental Science, Northwest Normal University. Water in the soil was 
extracted by an automatic vacuum condensation extraction system (LI-2200, BeiJingJianLing, 
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Beijing, China). The vacuum degree was controlled below 1 Pa/s, the heating temperature was 
105°C, and the extraction time was 3 h. After extraction, we randomly selected some samples and 
weighed them to ensure that the extraction efficiency was greater than 98%. The isotopic 
composition of the extracted soil water and precipitation was determined by a liquid water stable 
isotope analyzer (T-LWIA-45-EP, Los Gatos Research, California, USA). The measured 5!°O and 
5-H are expressed in thousandths, which are the difference relative to the Vienna Standard Mean 
Ocean Water. The equation can be expressed as follows: 
8 -[ 2-1} 0%, (1) 
standard 
where 6 is the isotopic composition (87H or 8!8O) in the water sample; Rsampie is the ratio of the 
heavy isotope to the light isotope in the water sample, i.e., 7H/'H and '80/!°O; and Rstandara is the 
isotope ratio in the Vienna Standard Mean Ocean Water. The standard error of measurement is 
+1.0%o for ôH and +0.3%o for 6!8O. Organic substances, such as methanol, ethanol, and other 
biovolatile compounds, may be mixed in the process of extracting soil water by low-temperature 
vacuum distillation, which will lead to large isotope differences in the process of laser analysis. 
To eliminate the pollution caused by organic substances, we used the spectral analysis software 
(Version 3.1.0.9, Los Gatos Research, California, USA) to establish the correction curve, and 
corrected the isotopes of soil water data based on the index of the measured absorption spectrum. 
The soil water content was calculated by the drying and weighing method. The wet soil was 
weighed with an electronic balance, baked for approximately 24 h in a constant temperature oven 
(105°C+2°C), dried to constant weight, and then weighed. The calculation formula is as follows: 


W, -W, 


2 0 


SWC = } 100% , (2) 
where SWC is the soil water content (%); Wi is the weight of wet soil in the aluminum box before 
drying (g); W2 is the weight of the aluminum box with dry soil after drying (g); and Wo is the 
weight of the empty aluminum box (g). 


2.3 Meteorological data 


The meteorological data during the sampling period were recorded using the automatic weather 
station in the meteorological field of Northwest Normal University (WatchDog 2000, Spectrum 
Technologies, Chicago, USA), including air temperature and relative humidity. Precipitation 
stable isotopes data are derived from precipitation collected in the meteorological field from 2019 
to 2021. The rain gauge was located next to the automatic weather station and was used to collect 
rainfall generated by each rainfall event from 2019 to 2021. 


2.4 Determination of the soil evaporation loss rate 


2.4.1 Line-conditioned excess (Ic-excess) coupled Rayleigh fractionation model 


In this study, we assumed that the fractionation of hydrogen and oxygen stable isotopes in the 
process of soil water evaporation conforms to the Rayleigh fractionation model (Clark and Fritz, 
2013). The soil evaporation loss rate (f; %) can be calculated as follows: 


R, =R,0- f)”, (3) 
where R, and R, represent the ratio of 7H/'H (or '80/'°O) to soil water and soil water source, 


respectively; and a represents the fractionation factor. 
Replace R, and R, in Equation 3 with òs and õp, respectively, we can obtain Equation 4: 


5, = (6, +1000) f)"” —1000, (4) 


where òs (%o) and 5, (%o) are the stable isotopes values of soil water and soil water source, 
respectively. 


a=(l-«,x10°)/a* , (5) 
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where a* is the equilibrium fractionation factor and & is the dynamic fractionation factor (Horita 
et al., 2008). 

Landwehr and Coplen (2006) proposed Ic-excess value to evaluate whether there is an isotopes 
shift between soil water and soil water source. In this study, lc-excess value represents the 
deviation degree of the isotopes value concerning the input water line. The equation of lc-excess 
value can be expressed as: 


Ic-excess=5°H — (ax 5'°O +b), (6) 


where a and b are the slope and intercept of the input water line, respectively. 

The physical meaning of Ic-excess is expressed as the degree of deviation of the isotopes value 
in the sample from the input water line, indicating the non-equilibrium dynamic fractionation 
process caused by evaporation. Therefore, Formula 7 can be obtained: 


le-excess = (°H, +1000)((1— f) )-1000-a( 0, +1000)(q— fy"*%?)-1000)-, (7) 


where 67H, and 8/80, represent stable isotopes of soil water. 
2.4.2 Craig-Gordon model 


Based on the principle of isotopes mass balance, we used the Craig-Gordon model to calculate the 
soil evaporation loss rate in an open liquid-vapor isotopes system (Gonfiantini, 1986; Sprenger et 
al., 2017). The formula of soil evaporation loss rate (f, %) is as follows: 


s-i- S82) , (8) 


where 5° is the limiting factor for isotopes enrichment (Gonfiantini, 1986). The formula is as 
follows (Gibson et al., 2014): 
+ + 
š = a B , (9) 
h-10 ` (e+e /a’) 
MN h-10°(é, +e" /a*) 
1-h+107e, 

where A is the relative humidity in the air; 5, is the isotopes value in the surrounding air; ¢* is the 
equilibrium fractionation coefficient (Gibson et al., 2008); and m is the enrichment slope. 54 can 
be calculated as follows: 


(10) 


T (11) 


We calculated a* and e* according to Horita et al. (1994): 


3 2 9 
ioil ct, e 2882 1620.17 TABAT 161,04 + 2999210 . 12) 
CH) 10 10 10° 
3 6 9 
10°In| eye, }=—7.685 + 2123810 1.666410  0:3504%10 | (13) 
(0O) T T 
e*=(a* -1)x1000, (14) 


where T is the temperature (K). & is the dynamic fractionation factor, and the formulae for 67H 
and ô!ŝO are as follows (Merlivat, 1978): 


Eem =X -A)x (1-0.9755), (15) 


Erdo) =nx(1—h) x (1—0.9723) , (16) 


where n is the aerodynamic parameter of the liquid-air interface during evaporation. Generally, it 
is 0.50 for saturated soil and 1.00 for dry soil (Benettin et al., 2018). Affected by precipitation 
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infiltration and soil evaporation, the topsoil is often in the process of alternating dry and wet 
conditions. Therefore, we set 0.75 as the average value of parameter n in this study, representing 
the average dry-wet condition of the soil (Mahindawansha et al., 2020). 

The slope of the local evaporation line (Sszz) is calculated as (Gonfiantini, 1986): 


[He ~§,)-(1+8, x0.001\(e, +e" a 
2 


1- h+ £, x 0.001 
SSEL = > (17) 
+ + 
Peoi O Ma e la J 
18 


1-h+e, x0.001 


where do is the isotopes value of the multi-year average precipitation. 
The intercept (IseL) of the soil evaporation line was calculated based on the isotopes of soil 
water (57H, and 8!80,): 


Iser = 5°H, —Ssep X 50, i (18) 
The isotopes of the soil water source (8!8Ointersect and 67Hintersect) is the intersection of the soil 
evaporation line and the input water line (Javaux et al., 2016). The formula is as follows: 


Oraa =. 


intersect 


(19) 


-5 SEL 


5H = að!$O +b. (20) 


intersect intersect 


2.5 Statistical analysis 


For a given mean value and standard error, we used Student's ¢ test to analyze the difference 
between the two mean values with unequal variances: 


paat Y (21) 
Jox toy’ 
where x and y are the two mean values; and ox and oy are their standard errors, respectively. The 
statistically significant difference was indicated by t>1.96 (P<0.05). 


3 Results 


3.1 Soil water content 


In terms of temporal variation, the soil water content of the plots without gravel-sand mulching 
changed significantly between different months (Fig. 2). The variation range was between 
0.17%/month and 3.00%/month. The average soil water content was higher in May and September, 
and relatively low in October, which was closely related to the water input during the sampling 
interval. The correlation coefficient between the soil water content and input water content was 
0.75 (P<0.01), with a positive correlation (Table 1). The monthly variation range of the soil water 
content in the plots with gravel-sand mulching was 0.86%/month to 2.94%/month, which was 
smaller than that in the plots without gravel-sand mulching. The maximum and minimum values of 
the average soil water content appeared in September and October, respectively. In terms of spatial 
variation, the soil water content at the soil depth of 0-30 cm was the most significant (P>0.05), 
that at the soil depth of 30-150 cm was relatively stable, and that at the soil depth of 100-150 cm 
was low. In short, the soil water content of the plots with gravel-sand mulching was higher than 
that of the plots without gravel-sand mulching both in time and space, indicating that gravel-sand 
mulching can indeed play a role in soil and water conservation. 


3.2 Stable isotopes of soil water 


During the whole sampling period, the average 5°H and 8!8O in the soil water source were 
—84.99%o and —12.21%o, respectively; the average 57H and 5!°O in the plots without gravel-sand 
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Fig. 2 Temporal and spatial variation of the soil water content. (a, c, and e), the plots with gravel-sand mulching; 
(b, d, and f), the plots without gravel-sand mulching. 


Table 1 Correlation analysis between the soil evaporation loss rate and other factors 


Factor feo) few Ficer RH T IW SWC Le-excess CW 
Fo 1.00 
few 0.97" 1.00 
Jac-excess) 0.98™ 0.98*™ 1.00 
RH —0.71* —0.66" 0.76" 1.00 
T 0.55" 0.50" 0.61" -0.92™ 1.00 
IW -0.44" -0.54" -0.54" - - 1.00 
SWC -0.73"" —0.70" -0.76" 0.29" —0.27" 0.75" 1.00 
Le-excess 0.56" —0.48"" —0.62™" 0.50"" -0.48"" 0.61” 0.87” 1.00 
CW 0.60" 0.77" 0.72" - - - - - 1.00 


Note: f('%o), the soil evaporation loss rate calculated by 6'°O using Craig-Gordon model; ffm), the soil evaporation loss rate 
calculated by ôH using Craig-Gordon model; fic-excess), the soil evaporation loss rate calculated by line-conditioned excess 
coupled Rayleigh fractionation model; RH, relative humidity; T, average temperature; IW, input water content; SWC, soil water 
content; Lc-excess, line-conditioned excess value; CW, crown width. *, P<0.05; “*, P<0.01; -, no data. 


mulching were —70.47%o and —9.46%o, respectively; and the average 8°H and 8!8O in the plots 
with gravel-sand mulching were —75.13%o0 and —10.49%o, respectively (Fig. 3). Since the water 
input of these plots are derived from more than the atmospheric precipitation, there will be 
irrigation in the early stage of plant growth, and we regard the irrigation water as strong 
"precipitation" and jointly fit the water line equation, which is called the input water line. 
Considering the input of precipitation and irrigation water in the later stage, we fit the input water 
line of each month and calculated the soil evaporation line of each month. The stable isotopes 
values of soil water in each month were lower than that of the input water line (Fig. 4), and the 
slope and intercept of the soil water evaporation line were also smaller than those of the input 
water line, indicating that there was no other water input in these plots and that the soil water was 
affected by evaporation. Regardless of whether the plots are mulched with gravel and sand, the 
stable isotopes values of soil water first decreased and then stabilized. The stable isotopes values 
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Fig. 3 Temporal and spatial changes of the stable isotopes values of soil water. (al-f1) represent ôH and (a2-f2) 
represent 5!8O. (a, c, and e), the plots with gravel-sand mulching; (b, d, and f), the plots without gravel-sand 


mulching. 
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of soil water at the soil depth of 0-30 cm changed most significantly (P<0.05) (the plots 
without gravel-sand mulching: 57H, —83.47%o— —30.27%o0; 5!8O, —11.96%o— —2.32%o; the plots 
with gravel-sand mulching: 57H, —85.13%0— —53.34%o; 5!8O, —12.10%o— —6.34%o). Moreover, 
the stable isotopes values of soil water at the soil depth of 0—30 cm were significantly richer 
than those at the soil depth of 30—150 cm, which is mainly due to the influence of evaporation. 
However, the stable isotopes enrichment of deep soil water in plot C and plot D was not 
affected by evaporation, which may be caused by the infiltration of enriched soil water in the 
surface layer. In terms of sampling time, there were some differences in the average stable 
isotopes values of each month, but there was no regular change with the sampling time. The 
monthly variation in the stable isotopes values of soil water in the plots without gravel-sand 
mulching (87H, 2.33%o/month—6.37%o/month; 5!8O, 0.15%o/month—0.87%o/month) was greater 
than that in the plots with gravel-sand mulching (87H, 1.89%o/month—3.89%o/month; 8!8O, 
0.19%0/month—0.59%o/month). The result was almost consistent with the monthly change of the 
soil water content, which indicated that gravel-sand mulching weakens the intensity of soil 
water evaporation. The depletion and enrichment of 57H and 5!8O were the results of 
evaporation, and the stable isotopes values in the plots with gravel-sand mulching were lower 
than those in the plots without gravel-sand mulching. Therefore, the change of the stable 
isotopes values of soil water is a good indication of soil evaporation under different mulching 
measures. 
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Fig. 4 Input water line, soil water evaporation line, and global meteoric water line in April (a), May (b), June (c), 
July (d), August (e), September (f), and October (g). The shaded areas are 95% confidence interval. 


3.3 Soil evaporation loss rate 


The evaporation fractionation signals of soil water stable isotopes at the soil depth of 0—30 cm 
were well preserved, so we discussed the monthly soil evaporation loss rate at this depth. We 
used the Ilc-excess value to express the deviation degree of the stable isotopes values of soil 
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water from the input water line; the smaller the value is, the stronger the evaporation 
fractionation is (Table 2). In terms of spatial variation, with the increase of soil depth, the 
Ic-excess values of the plots with and without gravel-sand mulching showed an increasing trend, 
indicating that the closer to the surface, the stronger the evaporation. In terms of temporal 
variation, the trend of Ic-excess value was nearly consistent with the change of the stable 
isotopes values of soil water. The average Ic-excess value of the plots with gravel-sand 
mulching was —4.86%o, and the average Ic-excess value of the plots without gravel-sand 
mulching was —7.79%o. Obviously, the gravel-sand mulching measure weakens the evaporation 
of soil water. 


Table 2 Line-conditioned excess (Ic-excess) values of the plots with and without gravel-sand mulching at the 
soil depth of 0-30 cm 


Pioi Soil depth Lc-excess (%o) 
(cm) Apr May Jun Jul Aug Sep Oct 
0-10 -14.68 -1.23 —10.87 -9.97 -8.11 -2.97 —4.84 
A 10-20 -7.43 —5.24 2.84 6.04 -5.13 -1.17 -3.55 
20-30 -5.74 —0.37 -3.67 -3.21 —4.04 —4.47 —4.85 
0-10 —16.34 —6.74 -7.33 -7.51 -5.54 -3.16 —4.64 
C 10-20 —6.95 -8.15 -2.01 -3.26 —0.70 -2.95 -5.74 
20-30 —4.73 —9.72 —0.87 -1.73 -1.57 —0.51 -5.27 
0-10 —18.28 —1.74 —10.82 -9.63 —6.29 —8.26 -9.76 
E 10-20 —9.90 -1.87 -5.16 -7.73 —4.15 —4.85 -3.39 
20-30 —4.86 -3.56 -1.77 6.94 -4.19 -3.13 —4.09 
0-10 -17.93 -1.61 -17.50 -13.39 -13.62 -5.75 -5.53 
B 10-20 —6.39 —9.73 -3.66 —6.67 -8.10 —4.30 -5.35 
20-30 -1.24 -3.34 -1.82 —4.99 -7.45 —4.46 —6.92 
0-10 -12.37 -8.11 -31.21 -12.51 —9.35 —6.21 -5.20 
D 10-20 -3.03 —8.94 -12.34 -5.50 -8.21 -3.65 —4.37 
20-30 -2.64 —5.94 —8.29 -3.55 —6.54 -2.33 —6.43 
0-10 —13.67 —0.63 —26.53 —18.53 -7.83 -11.99 —10.99 
F 10-20 -5.64 —0.01 -8.92 -9.43 —4.30 —8.20 -9.08 
20-30 —4.85 —0.10 —5.63 -9.50 -5.37 -3.39 -9.58 


Note: A, C, and E are the plots with gravel-sand mulching; B, D, and F are the plots without gravel-sand mulching. 


Figure 5 showed that, whether the soil is mulched with gravel and sand or not, evaporation 
was the strongest in April (the soil evaporation loss rate calculated by 57H using Craig-Gordon 
model (f7H)), 35.03%; the soil evaporation loss rate calculated by 5'8O using Craig-Gordon 
model (f('80)), 37.58%; and the soil evaporation loss rate calculated by the Ic-excess coupled 
Rayleigh fractionation model (fac-excess)), 58.31%) and relatively weak in May (fêm, 12.34%; 
feo), 22.07%; and fic-excess), 33.09%) and September (fêm, 10.73%; fdo), 16.37%; and fitc-excess)s 
30.83%). In the growing season of peach trees, the average soil evaporation loss rate in the 
plots with gravel-sand mulching was 11.54% for fêm, 17.19% for fdo), and 29.52% for 
Jic-excess), and the average soil evaporation loss rate in the plots without gravel-sand mulching 
was 21.69% (fêm), 27.96% (fto), and 43.94% (fac-excess)). Obviously, gravel-sand mulching 
reduces the consumption of ineffective water. The difference in the soil evaporation loss rate 
between the plots with and without gravel-sand mulching was the largest in August, which 
indicated that the effect of gravel-sand mulching measures is most prominent in August, so we 
can increase the thickness of gravel-sand mulching in this month. 
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Fig. 5 Temporal variation of the soil evaporation loss rate (a), temperature, relative humidity, soil evaporation 
line slope (Sszz), and precipitation (b). fH), the soil evaporation loss rate calculated by 57H using Craig-Gordon 
model; filo), the soil evaporation loss rate calculated by 5'8O using Craig-Gordon model; fic-excess), the soil 
evaporation loss rate calculated by line-conditioned excess coupled Rayleigh fractionation model; the subscripts 
of mulching and non-mulching represent the plots with and without gravel-sand mulching, respectively. The 
shaded area is the soil evaporation loss rate in similar environments, and the yellow dashed line is the average soil 
evaporation loss rate in similar environments (Schlesinger and Jasechko, 2014). 


4 Discussion 
4.1 Influencing factors of soil evaporation 


We applied the Ic-excess value to qualitatively determine the intensity of soil evaporation, and 
used the lIc-excess coupled Rayleigh fractionation model and Craig-Gordon model to 
quantitatively calculate the soil evaporation loss rate. However, soil evaporation is the result of 
many factors. The lc-excess value includes many other factors that affect soil evaporation when it 
expresses the intensity of evaporative fractionation, therefore, we analyze the correlation between 
Ic-excess value and other factors (Table 1). The Ic-excess value was positively correlated with 
relative humidity, with a correlation coefficient of 0.50, and negatively correlated with average 
temperature, with a correlation coefficient of —0.48; and the correlation coefficients of lc-excess 
value with input water content and the soil water content were 0.61 and 0.87, respectively. The 
Ic-excess value can reflect the strength of evaporation. The results showed that evaporation is 
stronger in April, which may be because the branches and leaves of the peach trees are not fully 
developed at the early stage of growth. Therefore, we performed a correlation analysis between 
crown width and evaporation, and the average correlation coefficient was 0.70 (P<0.01; Table 1). 
In other words, lusher branches and leaves can reduce the evaporation of soil water to a certain 
extent. The relatively weak evaporation in May and September is not affected by crown width, 
which may be attributed to the amount of water input (Hasselquist et al., 2018; Dubbert et al., 
2019). There is a negative correlation between evaporation and input water content, and the 
average correlation coefficient is —0.51 (P<0.05; Table 1). Due to the existence of irrigation in 
early May, the total water input in May is larger than that in September, so the soil evaporation 
loss rate in May is slightly smaller than that in September. The average correlation coefficients of 
evaporation with relative humidity, average temperature, and soil water content were —0.71, 0.55, 
and —0.73, respectively (P<0.01; Table 1). The soil evaporation loss rate is negatively correlated 
with relative humidity and positively correlated with average temperature, which is consistent 
with the fact that evaporation will be weaker with higher humidity and stronger with higher 
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temperature. These findings are consistent with the research results of Xiang et al. (2021), which 
showed that soil evaporation is affected by multiple factors. 


4.2 Calculation of the soil evaporation loss rate 


Soil evaporation will cause water isotopes fractionation, resulting in changes of soil water stable 
isotopes values (Zimmermann et al., 1966; Allison and Barnes 1983). Therefore, the influence of 
soil evaporation can be explored through the changes of soil water stable isotopes values 
(Sprenger et al., 2017). However, it is generally believed that the direct evaporation depth can 
occur only at the top of soil layer (Or et al., 2013; Piri et al., 2020), but some studies showed that 
the fractionation signal of isotopes can reach the soil depth of 20—30 cm in temperate regions, and 
it can reach the soil depth of 50 cm in Mediterranean climate regions (Sprenger et al., 2016). In 
this study, the fractionation signals of soil water isotopes at the soil depth of 0-30 cm are 
prominent, so we calculated the soil evaporation loss rate at this depth. 

There are two methods to calculate the proportion of water entering the soil returning to the 
atmosphere through soil evaporation (i.e., the soil evaporative loss rate). The first method is based 
on the Rayleigh fractionation (Clark and Fritz, 1997), and the second method is based on the 
Craig-Gordon model (Gonfiantini., 1986; Sprenger et al., 2017). To highlight the isotopic 
advantage, we used the Ic-excess coupled Rayleigh fractionation model in this study. In general, 
plot A and plot B in the Ic-excess coupled Rayleigh fractionation model represent the slope and 
intercept of the local precipitation line, respectively. Since there is irrigation in these plots, we 
fitted the input water line of each month, and the Ic-excess value represents the deviation degree 
of isotopes value in the sample from the input water line, which makes the Ic-excess coupled 
Rayleigh fractionation model more accurate. When the Craig-Gordon model is applied to 
calculate the soil evaporation loss rate, the soil water source needs to be determined, which is 
usually determined by the soil evaporation line and local meteoric water line, and the intersection 
of the two is the soil water source (Javaux et al., 2016). Most studies only identified one soil 
water source to calculate the soil evaporation loss rate in time series (Gibson et al., 2016; Che et 
al., 2020; Mahindawansha et al., 2020), ignoring the input of precipitation and irrigation water. 
For this reason, we determined the soil water source of each month through the input water line 
and soil evaporation line of each month (Fig. 4 and Table 3), which can improve the 
computational accuracy of the Craig-Gordon model to a certain extent. 


Table 3 Soil water source during April—October in 2021 


Soil water source 


Month 


53O (%o) ôH (%o) 
April 2021 —11.45 -79.86 
May 2021 -12.50 -87.47 
June 2021 -12.35 —85.99 
July 2021 -11.79 —81.65 
August 2021 -12.81 —89.20 
September 2021 -12.20 —84.67 
October 2021 -12.39 -85.92 


We calculated the soil evaporation loss rate at the soil depth of 0-30 cm by the Ic-excess 
coupled Rayleigh fractionation model and Craig-Gordon model. The change trends of fic-excess), 
feo), and fêm with time are consistent, and the correlation coefficient of the lc-excess value with 
fac-excess), fto), and fêm are —0.62, —0.56, and —0.48, respectively. Then, we compared the soil 
evaporation loss rate in similar environments (Fig. 5), and found that fac-excess) is completely 
located in the shaded area, followed by f('8o), which is consistent with the correlation between 
evaporation and the Ic-excess value. The relatively poor calculation result of fêm) may be because 
°H is more affected by artificial intervention such as fertilization in cultivated soil, while 18O is 
less affected by artificial intervention. Moreover, when the Ic-excess value was used to calculate 
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evaporation, the influence of 7H and '8O was weakened, so fac-excess) achieved the best result, 
followed by f('8o). 


4.3 Effects of gravel-sand mulching on soil evaporation 


We compared the difference of the soil evaporation loss rate at the soil depth of 0-30 cm. 
Throughout the growing season, the evaporation of soil water in the plots with gravel-sand 
mulching was significantly lower than that in the plots without gravel-sand mulching. In addition, 
the soil water content in the plots with gravel-sand mulching was also higher than that in the plots 
without gravel-sand mulching. The difference of the soil evaporation loss rate caused by 
gravel-sand mulching is the largest in August; that is, the effect of gravel-sand mulching on 
inhibiting soil evaporation is most prominent in August. 

Gravel-sand mulching is the main reason for the difference of soil water evaporation signals. 
Surface mulching changes the amount of water and heat exchange between the Earth and 
atmosphere, affecting the soil water phase transition (Fu et al., 2015), soil water and heat transfer 
(Lu et al., 2019), and soil water evaporation (Chen et al., 2019). In agriculture, gravel-sand 
mulching is often used to reduce the ineffective evaporation of soil water during crop growth 
season (Tan et al., 2019). This approach is used because the surface mulched with gravel and sand 
can cut off the rising path of the lower soil moisture and change the position of the evaporation 
interface, so that the rising soil moisture is stored near the surface and any water vapor must pass 
through the sand layer to escape into the atmosphere. Gravel-sand mulching avoids partial soil 
water loss due to evaporation and reduces the soil evaporation loss rate. The study of farmland 
water status under the condition of gravel-sand mulching will be helpful to guide production 
practices and further improve the basic theory of mulching technology in agricultural water-saving 
to alleviate the substantial pressure caused by the shortage of water resources in Northwest China, 
which has important scientific value and practical significance for guiding the efficient utilization 
of agricultural water resources. Therefore, we suggest to implement gravel-sand mulching 
measures in the crops growing season in water-deficient areas, and appropriately increase the 
thickness of gravel-sand mulching in July and August to achieve better water retention. 


5 Conclusions 


We monitored the variation characteristics of soil water content and stable isotopes values during 
the growth season of Baifeng peach (Amygdalus persica L.) in 2021 and analyzed the 
characteristics of soil water and the soil evaporative loss rate when mulched with and without 
gravel and sand. Our results showed that the average soil water content in the plots with 
gravel-sand mulching is 1.86% higher than that without gravel-sand mulching, and the monthly 
variation of the soil water content is smaller in the plots with gravel-sand mulching than that in 
the plots without gravel-sand mulching. Second, the Ic-excess value was negative for both the 
plots with and without gravel-sand mulching. The average Ic-excess value of the plots without 
gravel-sand mulching is 2.93% lower than that of the plots with gravel-sand mulching. The 
Ic-excess value has good correlation with relative humidity, average temperature, input water 
content, and soil water content. The difference of the soil evaporation loss rate between the plots 
with and without gravel-sand mulching is the largest in August, and the inhibition effect of 
gravel-sand mulching on soil evaporation is most prominent in August. In conclusion, we 
quantified the influence of gravel-sand mulching on the soil water content and soil evaporation 
loss rate in this study, which provides some basic data support for improving water use efficiency 
and saving water resources. 
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